Prom laser cooling to aging: a unified Levy flight description 
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Intriguing phenomena such as subrecoil laser cooling of atoms, or aging phenomenon in glasses, 
have in common that the systems considered do not reach a steady-state during the experiments, 
although the experimental time scales are very large compared to the microscopic ones. We revisit 
some standard models describing these phenomena, and reformulate them in a unified framework 
in terms of lifetimes of the microscopic states of the system. A universal dynamical mechanism 
emerges, leading to a generic time-dependent distribution of lifetimes, independently of the physical 
situation considered. 



I. INTRODUCTION 

One of the important paradigmatic changes in the 
twentieth century physics was the introduction of proba- 
bilities, and hence fluctuations, as a central concept in the 
description of the physical world. This arose through the 
development of both quantum mechanics and statistical 
mechanics. Along this line of thought, a physical observ- 
able is often decomposed into an average value plus some 
fluctuations. Accordingly, fluctuations are most often 
seen as small perturbations around a well-defined mean 
value, meaning that the probability distribution of the 
observable is narrow. 

However, it has become clear since the mathematical 
work of Paul Levy in the 1930'Sfi^ that interesting sit- 
uations arise when the probability distribution becomes 
so broad that the average value of the observable is not 
even defined (i.e., it is formally infinite). This appears 
in particular when the tail of the distribution p(x) of the 
variable x behaves as a power law 



P(x) 



,1+a ' 



00, 



(1) 



with < a < 1. In the early 1970's, it was realized 
that Levy's results could be relevant to the study of ran- 
dom walksj 3 -^ with physical applications for instance in 
disordered conductors^ If the broadly distributed vari- 
able corresponds to the sojourn times in some micro- 
scopic states of a system, some interesting and exper- 
imentally testable properties arise, like anomalous dif- 
fusion, or aging. In the last fifteen years, this topic 
has become more and more popular in the physics re- 
search community. Broad distributions have found ap- 
plications in physical processes as diverse as turbulent 
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ran- 



and chaotic transport, PjLs. glassy dynamics 
dom walks in solutions of micelles^ diffusion of spec- 
tral lines in disordered solidS ) 13 i 14 fluorescence of single 
nanocrystalsji^ laser trapped ions ; 16 ' 17 or laser cooling of 
atomic gase o 18 i 19 (for which the Nobel Prize was awarded 
in 1997), to quote only some of them. 

The aim of this paper is to present a brief overview, 
within a unified formalism, of some elementary physical 
models in which the dynamics involves broad distribu- 
tions of sojourn times, and to highlight the surprising 



phenomena that may appear from such a simple dynam- 
ics. Using simple probabilistic arguments, we show that 
the equilibrium distribution becomes non normalizablc 
in a given parameter range, thus ruling out the existence 
of an equilibrium state. We then describe the dynam- 
ics of the different models in terms of lifetime variables 
-the mean time spent in a given microscopic state. This 
leads to a universal characterization of systems with non- 
stationary dynamics at large time. In this description, 
the common feature of these systems is the presence of 
power-law tails in the distribution of lifetimes of the mi- 
croscopic states, leading to an infinite average lifetime. 
We also compute the dynamical probability of lifetimes, 
that is, the probability for the system to occupy at time 
t a state with a given lifetime, and find at large time a 
universal form for this distribution. 

The paper is organized as follows. In Sect. HU we give 
two examples of physical situations where the time spent 
by the system in different microscopic states becomes 
broadly distributed, and present simple models aimed at 
describing these physical situations. In Sect. IIIH we show 
that the two models can be described within a unified 
framework, and compute the corresponding dynamical 
distribution of probability, before discussing some phys- 
ical phenomena that can be explained by these results. 
Finally, some technical details are reported in Appen- 
dices [3] and [BJ and we suggest in Appendix [Cl a project 
appropriate for students. 



II. SIMPLE PHYSICAL EXAMPLES 

A. Laser cooling of atoms 

Laser cooling of atomic gases consists in reducing the 
momentum spread of atoms thanks to momentum ex- 
changes between atoms and photons. The so-called "sub- 
recoil" laser cooling consists in reducing the momentum 
spread of the atoms to less than a single photon momen- 
tum hk, where k = 2t:/X is the wavevector of the light, 
and h is the reduced Planck constant. This cooling is 
achieved by introducing a momentum dependence in the 
photon scattering rate, in such a way that it decreases 
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strongly in the vicinity of p = 0, where p is the magni- 
tude of the atomic momentum p. When an atom reaches 
by chance this region of momentum space, it tends to 
stay there for a long time, since the photon scattering 
rate is very low. This leads to an accumulation of atoms 
at small momenta, that is, to a cooling of the atoms i 18 ^ 19 
The mechanism of subrecoil laser cooling is illustrated 
in Fig. [TJ Anytime a photon is absorbed and sponta- 
neously reemitted by an atom, the atomic momentum 
undergoes a random jump of the order of hk, since spon- 
taneous emission occurs in a random direction. Thus, 
the repetition of absorption-spontaneous emission cycles 
generates a random walk of the momentum of the atom, 
with momentum-dependent time intervals between two 
jumps. Following this random motion, the atom may 
end up in the low momentum region, where it stays for 
a long time, leading to the accumulation of atoms men- 
tioned above. 

The time interval r s between two jumps is called the 
sojourn time. It depends on the momentum of the atoms, 
and it also fluctuates from one realization of the process 
to the other, since photon scatterings occur in a random 
way. An important notion to characterize the dynamics 
is the mean sojourn time at momentum p, that is the av- 
erage time spent at momentum p between two successive 
photon scatterings. This mean sojourn time is denoted 
as the lifetime r(p), and it takes a well-defined value for 
each momentum p. The lifetime r(p) is the inverse of the 
photon scattering rate R(p), so that the determination 
of the lifetime requires some knowledge about the scat- 
tering rate. Subrecoil laser cooling methods are based 
on atom-light interaction processes such that the photon 
scattering rate R(p) vanishes for p = 0. In most cases, 
R(p) is a regular function of p 2 — p 2 due to rotational 
invariance, so that the small p behavior of R(p) is given 
by 

R(p)=r^(^)\ p^O, (2) 

\Pm J 

where tq is a time scale, and p m is a momentum scale, to 
be specified below. The exponent (3 is an even integer, 
(3 = 2, 4, or 6, which depends on the specific physical 
mechanisms used in the experiment j 1 8 1 1 9 It then follows 
that r(p) behaves at small p as a (diverging) power law, 

r(p) = ro (y)^ p ^°- (3) 

Hence for small enough p, the time spent at momentum 
p becomes very large with respect to the time spent at 
large momentum. To be more specific, one can find a mo- 
mentum scale p m <C hk (where hk is the typical change 
in atomic momentum when a photon is scattered) such 
that R(jp) behaves as a power law for p < p m , and such 
that the time spent by the atom in the region p > p m 
is very small as compared to the time spent at p < p m . 
Note that, as p m <C hk, values of p smaller than p m 




FIG. 1: Schematic picture of the subrecoil laser cooling pro- 
cess in momentum space. Photon scatterings induce a random 
walk for the momentum p of the atoms (p is the magnitude 
of p ), with steps of the order of hk, the photon momentum. 
The mean sojourn time r(p) at momentum p becomes very- 
large for small momenta (i.e. for p <C Pm), thus generating an 
accumulation of atoms in the vicinity of p = 0. 

are reached in an essentially uniform way, after a photon 
scattering from a larger momentum p w hk. 

Accordingly, to model in a simplified way subrecoil 
laser cooling, one neglects the time spent outside the re- 
gion p < p m , and simply assumes that the atom leaving 
this low-momentum region instantaneously comes back 
to it, choosing at random a new momentum within the 
sphere of radius p m , with a uniform distribution in d di- 
mensions. Thus the distribution p(p) of the momentum 
p reached after a jump into the low-momentum region is 

p(p) = ^Lp d -\ 0<p< Pm . (4) 

Pm 

Note that the absence of memory between two successive 
values of the momentum means that the dynamics chosen 
here corresponds to a Markov process -see Sect. IIII Al 

Let us now introduce another probability distribution, 
namely the probability P s t (p) for an atom to have a mo- 
mentum p at a given time, if a steady state is reached. 
This probability is proportional to the lifetime r(p), 
times the probability p(p) to reach a momentum p af- 
ter a jump. Hence, P s t{p) reads 

P st (p) = jL- p(p)r(p) = T -fkp-Pp{p), (5) 

where the average lifetime {t) p , computed over the dis- 
tribution p(p), is given by 

(r) p = dp p(p) T (p) = dT 0P t d dpp^ 1 ^. (6) 
Jo Jo 

The integral in Eq. © converges only if d > f3. For 
d < (3, the average lifetime is infinite, so that the distri- 
bution P s t(p) given in Eq. ([5]) is not normalizable. Thus 
no stationary state can be reached, and the system con- 
tinues to cool for infinitely long times. Accordingly, the 
divergence of the average lifetime turns out to play a key 
role in the phenomenon of subrecoil laser cooling, as we 
shall see in more details in Sect. IIII1 Note also that the 
divergence of (t) p explains why direct treatments of sub- 
recoil laser cooling using the equations of quantum optics 
become very difficult 
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FIG. 2: Schematic view of the energy landscape in the trap 
model: a particle is trapped within deep energy wells, or 
traps, of energy U — —E (E is called the depth of the trap) 
and can only escape through thermal activation. Once the 
particle leaves the current trap, it chooses another one at 
random among all the others. 



B. Aging and trap dynamics 

Glasses are materials for which the relaxation time 
becomes, at low temperature, much larger than the 
accessible experimental time scales. The phenomenol- 
ogy of glasses has been shown to be very rich, includ- 
ing in particular non-exponential relaxations and aging 
phenomena ; 20 ' 21 in which case the relaxation time of the 
system is proportional to its age (i.e., the time spent in 
the low temperature phase). 

A simple model for glassy dynamics, that gives inter- 
esting insights into the underlying physical mechanisms, 
is the trap model. In this model, the system is repre- 
sented by a particle in a high dimensional phase space 
and is trapped into low energy states (traps), as schemat- 
ically shown in Fig. [51 These traps are characterized by 
their depth E (i.e., the energy barrier to be overcome so 
as to leave the trap), considered to be random variables 
with distribution p(E). The particle can leave a trap by 
thermal activation only, and once a trap is left, a new one 
is chosen at random, with uniform probability. In physi- 
cal terms, this means that the time spent by the system 
in high energy states, that are not explicitely described 
in the model, is very short. In the simplest case, the dis- 
tribution p(E) is taken to be exponential (which can be 
justified on the basis of extreme value statistics 2 ^), 



P{E) 



k B T q 



,-E/k B T g 



E>0. 



(7) 



where k B is the Boltzmann constant. Up to now, k B T g 
is simply an energy scale. In the following, we consider 
the limiting case of an infinite number of traps. 

The dynamics of the model is defined by giving the dis- 
tribution ip(r s \E) of the sojourn times r s in a trap with 
depth E. We consider here a Markov process, which 
means that sojourn times associated to successively vis- 
ited traps are uncorrelated. For Markov processes, the 
distribution of sojourn times is exponential, as discussed 
in Sect. IIII A1 We thus consider the following exponential 



form for ip(r s \E), namely 
<p(r.\E) = - 



-(E) 



,-t s /t(E) 



(8) 



where t(E) is the mean sojourn time, i.e., the trap life- 
time. Note that the probability per unit time to es- 
cape the trap is constant in a given trap, and is equal 
to 1/t(E). 

From thermal activation arguments, the trap lifetime 
t(E), is given by the usual Arrhenius law, 



-(E) = t e 



(9) 



To being a microscopic time scale that characterizes the 
random motion of the particle within the trap. 

In steady state, the probability to be in a trap of depth 
E is proportional to the trap lifetime, times the density 
distribution p(E) of traps of depth E, leading to 

P eq (E) = p(E) r(E) = p(E) e E ' k * T (10) 



{r)i 



which is nothing but the usual Boltzmann-Gibbs distri- 
bution (note the + sign in the exponential, as E is the 
opposite of the energy U of the trap) . The average life- 
time (t)e is given by 



/>oo 

(t) b =t dEp(E)e E ^ T , 
Jo 



that is, 



(t)e = 



TO 



k B T„ 



POO 








/ c\E exp 
'o 




~T gJ 


— 

k B _ 



(11) 



(12) 



If T < T g , the equilibrium distribution P eq (E) ceases to 
be normalizable since the integral defining (r) B diverges. 
Thus here again, the divergence of the average lifetime 
turns out to play a major role, ruling out the existence 
of a steady state and leading to a non-stationary dynam- 
ics at large times. Accordingly, starting from a uniform 
distribution among traps (to be interpreted as an infi- 
nite temperature initial condition), the average dynami- 
cal energy of the system slowly drifts toward lower values, 
without ever being able to reach a steady state. This will 
appear explicitly in the next section. 



III. UNIFIED DESCRIPTION IN TERMS OF 
LIFETIMES 

A. Distribution of lifetimes 

We have seen that broad distributions of lifetimes can 
be found in different physical problems belonging to 
different fields of physics. These problems are usually 
studied using the natural physical variable, like momen- 
tum for laser cooling^ or energy barrier for the trap 
modeL 23 ' 24 Yet, once expressed in terms of lifetimes, 
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these problems become equivalent from a formal point of 
view, at least when considering the dynamical distribu- 
tion of lifetimes. Differences only appear when one wants 
to look at the predictions for the physical variables, as 
discussed in Sect. MIDI 

Let us compute the distributions of lifetimes for 
the laser cooling model as well as for the trap model. 
Starting from Eqs. ([3]) and (U|) -resp. from Eqs. ([7]) 
and ([9])- and using the relation p(p)\dp\ = ip(r)\d.T\ - 
resp. p(E)\dE\ — ^(T)\dr\- one finds in both models a 
power law behavior for ip{r), 

A 

VK T ) = ^T+^, T > Tmin, (13) 

with a = d/P for laser cooling, and a = T/T g for the trap 
model; A and T m ; n are model-dependent parameters, the 
precise value of which being unimportant to the following 
discussion. For a < 1, the first moment (r) diverges, 
which means that the lifetime averaged over all possible 
states is formally infinite. This has deep consequences on 
the dynamics, as discussed below, since the system is in 
this case no longer able to reach a steady state. 

At this stage, let us emphasize that two different no- 
tions of time intervals appear in these problems, namely 
the sojourn time t s and the lifetime r. The sojourn 
time is the time spent between two successive events - 
scatterings by a photon for laser cooling, or escapes from 
a trap in glassy dynamics. The lifetime is associated to 
a microscopic 'state' of the system, and is defined as the 
average value of the sojourn time in this state. Lifetimes 
can also be used to label the corresponding states, instead 
of using the natural physical variables of the problem. In 
addition, it is important to notice that the models con- 
sidered here are Markov processes. This means that for a 
given lifetime r, the probability to exit the current state 
during a short time interval dt is dt/r. As a result, the 
distribution lp{t s \t) of sojourn times t s in a state with 
lifetime r is exponential for Markov processes, 

cp(r s \T) = -e-^ T . (14) 

T 

Then the full distribution ?/> s (t s ) of sojourn times, aver- 
aged over all microscopic states, is 

f°° Z" 00 dT 

Mt.)= <1t<p(t s \t)iP(t)= —^(j) e -^l\ 
Jo Jo T 

(15) 

Note that if ip(r) has a power law tail for r — > oo, then 
ip s (t s ) also has a power law tail, with the same exponent 
as tp(r). 

B. Master equation for the dynamical distribution 
of lifetimes 

To obtain a unified picture for the physical models pre- 
sented in Sect.|TTl we characterize the states of the system 
by their lifetimes r, without any reference to the physical 



meaning of these states. We wish to study the probabil- 
ity density P(t, t) to be in a state with lifetime r at time 
t. For a Markov process, the evolution of P(r, t) is given 
by the following master equation (see Appendix [A")) : 

dP f°° 

— (r,i) = -P(r,t) J dr'W(T^r') (16) 

p OO 

+ / At' W{t' ->t)P(tV), 
Jo 

where W(t — > r') is the transition rate from states with 
lifetime r to states with lifetime r'. The master equation 
is simply a balance equation for the probability P(t, i), 
taking into account the probability to leave states with 
lifetime r, and the probability to reach them starting 
from another state. Consistently with the models pre- 
sented in Sect. HH we make the important assumption 
that the lifetimes r l5 t 2 , ... r n of the states succes- 
sively occupied by the system are statistically indepen- 
dent. Hence, at each transition, a new lifetime r' is cho- 
sen at random according to the distribution ip( T ')- Given 
that the probability per unit time to leave state r is 1/r, 
the transition rate thus takes the form 

W(t^t') = -Mt'). (17) 
r 

The lack of correlation between the states before and 
after a jump occurs physically because of the strongly 
stochastic underlying dynamics. In subrecoil laser cool- 
ing, it comes from the fact that, in between two trap- 
ping events, the atoms actually perform a random walk 
outside the trap, which takes a negligible time, but un- 
correlates two consecutive trapping momenta. A simi- 
lar scenario holds for aging dynamics, where the system 
wanders among high energy states and fully decorrelates 
before finding a new trap. 

Note that the independence property, which is stronger 
than the Markovian property, would lead to a trivial 
stationary process if all states had essentially the same 
lifetime. However, as shown below, interesting non- 
stationary effects appear when the lifetimes are broadly 
distributed. 

Using the form of the transition rates given in Eq. (fT7|) 
and the normalization of ip( T ')i the master equation 
Eq. (flTjl) becomes 

dP 1 f°° dr' 

— (r,t) = --P(T,t)+^(T)j o —P(r',t). (18) 

C. Solution of the master equation 

The master equation Eq. (fT8]) can be solved exactly 
using a Laplace Transform with respect to time t. Cal- 
culations are described in details in Appendix [Bl and we 
summarize the main results in the following. The Laplace 
Transform P(r, s) of P(t, t) is defined as 

/>oo 

P(r,s) ee / dt e~ st P{T,t). (19) 
Jo 
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If we assume that at time t — 0, the system is distributed 
completely at random over all possible states, then the 
initial distribution P(r,t = 0) reduces to ip(r). In this 
case, P(t, s) can be written in a simple form 



P(r,s) 



1 



sr*(s) 1 + st' 



(20) 



where the time scale r*(s), to be interpreted as an effec- 
tive average lifetime, is defined as 



:(s) 



dr 



(21) 



Note that P(t,s) is normalized to 1/s with respect to 
t. It should be emphasized that Eq. ([2H)l is an exact 
result which does not rely on some long time regime (s <C 
Tq 1 ) approximation. The influence on the dynamics of 
the finiteness or divergence of the average lifetime (r) 
now appears clearly, since (r) = lim s _vo T l (s) ■ If ( r ) is 
finite, the distribution P(r, s) converges for s — > (to 
be understood physically as t — > oo) to the equilibrium 
distribution (up to the usual 1/s factor) 



P(t,s) 



a(r) 



rtp(r), 



(22) 



Thus in this case, the system reaches an equilibrium state 
at infinite time. Note however that the convergence to 
this asymptotic state may be very slow, typically as a 
power law with an exponent which becomes small when 
(r) is large. 

In contrast, if (r) is infinite, the effective average life- 
time r*(s) grows without bound when s — ► 0, and no 
stationary state is ever reached. If the lifetime distribu- 
tion ip(r) has a power law tail, 



1>(t) 



A+a ■ 



(23) 



with < a < 1, then the out-of-equilibrium solution of 
the master equation exhibits some interesting non trivial 
properties, in particular a scaling form at long times, that 
is for s — > 0. To be more specific, the distribution P(t, s) 
can be written asymptotically as 



P(T, S) « <f,(6 

where the scaling function 
Eq. [(20")). namely 



0. 



(24) 



(x) can be determined from 



1 



sin(7ra) x a (\ + x) 



(25) 



To come back to physical time, the inverse Laplace Trans- 
form of P(t, s) in the scaling regime can be computed. 
As expected, the resulting distribution P(t, t) also takes 
a scaling form, 



P(r,t) 



Ml)- 



(26) 



which shows that the typical momentum r t at time t 
is of the order of t itself. The scaling function 4>(u) is 
such that 4>(u) ~ u~ a for u — > and <p(u) ~ for 
u — > oo (note that </> is not the Laplace Transform of </>). 
The asymptotic behaviour of P(t, t) for t < ( or r > t 
can then be computed, and is found to be qualitatively 
the same as that of P(t, s = l/t). Neglecting multiplica- 
tive constants, the distribution P(t, t) takes for r -C t a 
pseudo-equilibrium form, 



P(r,t)~r^(r) J 



r < t, 



(27) 



whereas for r ^> it remains proportional to the a priori 
distribution V>( T )> 



P(r,t)^^(r), 



r > t. 



(28) 



This means that states with small enough lifetimes have 
been visited a large number of times and are essentially 
equilibrated. In contrast, states with very large lifetimes 
have been visited only once, and the precise value of their 
lifetime is not felt by the dynamics at this stage. 



D. Link between lifetime and physical observables 

Even though the different models introduced in Sect.lTTl 
share a common dynamical distribution P(r,t), physical 
observables do not behave in the same way in these mod- 
els. The reason for this is that the functional relations 
between these observables and the lifetimes differ much 
from one model to the other. 

In the subrecoil laser cooling problem, the relation be- 
tween momentum p and lifetime r is p — p m (T/ro) _1///3 . 
Hence from Eq. (|26|) . the dynamical distribution of mo- 
mentum P m (p, t) can be written, in the long time regime, 
as 



P m (p,t) 



tp+P \ tpP 



(29) 



Thus the distribution of momentum becomes more and 
more narrow, and the typical momentum p*(t) reached 
after a large time t is given by 



P*(t)~p m {t/T )- 1 " 3 . 



(30) 



which shows that the atoms can, in principle, be cooled 
down to arbitrarily low temperatures, provided that one 
waits long enough, and that no other effects are present. 
The above result is indeed in remarkable agreement with 
the observations made in laser cooling experiments, 19 
and leads to a clear understanding of this intriguing phe- 
nomenon, of purely dynamical origin. 

In the trap model of glass aging, the energy E is related 
to the lifetime r through E = — T\sx(t/tq). Since the 
typical lifetime r t at time t is of the order of t, the average 
energy E*(t) behaves as 



E*(t) w -Tln(iAo). 



(31) 
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Such a logarithmic behavior of the energy is indeed fre- 
quently observed for aging systems, and is accounted for 
in a simple way in the present approach. Indeed, the lat- 
ter gives a precise meaning to the word 'aging', as it is 
clear from the scaling form Eq. (|26p and from the average 
energy given in Eq. (j3Tj) that the system keeps track of 
its age. 

IV. CONCLUSION 

In this note, we illustrated on two simple examples 
how the apparently formal broad distributions can lead 
to deep insights and to a unified view on the physics of 
systems that do not reach an equilibrium state. This for- 
malism leads in a straightforward way to the long-time 
behavior of such systems, and shows that this behavior 
is universal, only depending on the exponent a of the 
power-law tail of the lifetime distribution. In addition, 
this approach gives a well-defined meaning to the notion 
of "quasi-cquilibrium" at large times, which is often re- 
ferred to in such systems. In a quasi-equilibrated system 
at a large time t, microscopic states with lifetimes less 
than t are essentially equilibrated, whereas those with 
a lifetime greater than t are not. Finally, let us men- 
tion that more general situations, for which the return 
time between two trapping events is not negligible, or 
even broadly distributed (while this time was neglected 
in the present approach), may also be considered. Such a 
generalization does not change the main qualitative con- 
clusions reported here^ 
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APPENDIX A: MASTER EQUATION 
FORMALISM 

On general grounds, let us consider a system described 
by a finite set of configurations (or microscopic states) 
{C}. The set {C} may be called the configuration space of 
the system. We assume that the system evolves according 
to a markovian stochastic dynamics, which means that 
the system may change configuration randomly in time, 
and without memory of the configurations occupied ear- 
lier. Such a random dynamics is often called a Markov 
process. To be more specific, let us assume that the sys- 
tem is in an arbitrary configuration C at time t. The 
probability that the system goes from C to a new config- 
uration C between time t and t + dt, where dt is an in- 
finitesimal time interval, is equal to W(C — » C')dt, which 



defines the transition rate (i.e., probability per unit time) 
from C to C, W(C — ► C). Roughly speaking, the marko- 
vian property of the dynamics states that W(C — > C) 
depends neither on the history of the system before ar- 
riving at configuration C, nor on the time spent at config- 
uration C. Let us also emphasize that W(C — > C) should 
be understood as the probability per unit time to go to 
configuration C given that the system is in configuration 
C; it is thus a conditional probability. 

Let us now introduce the probability Pc (t) for the sys- 
tem to be in a configuration C at time t. The evolution 
with time of Pc(t) is described by a master equation, 
which is nothing but a balance equation for the proba- 
bility. It may be obtained in the following way. In a very 
short time interval dt, the variation dPc of Pc(t) is due 
on the one hand, to the probability rf$g ut to leave the 
configuration C during the time interval dt, and on the 
other hand to the probability to reach the configu- 
ration C from any other configuration C in the same time 
interval. This yields 

dP c = -d$£ ut + d$j?. (Al) 

These two contributions can be evaluted as follows. 
Starting from C, the system may go to any other configu- 
ration C, so that the corresponding probability is simply 
the sum of all such possible transition probabilities. The 
probability to be in configuration C and to go to config- 
uration C within dt is W(C -> C')P c {t)dt, leading to 

d$° ut = W(C -» C')P c {t)dt. (A2) 

In a similar way, dQ™ is obtained by summing the transi- 
tion probabilities from all other configurations C to con- 
figuration C, 

d<f>£ = W(C -> C)P e >{t)dt. (A3) 
Substituting into Eq. (|A1|) . one obtains 

(A4) 

which is precisely the master equation describing the 
stochastic evolution of the system. Note that for simplic- 
ity, the above argument was restricted to systems with a 
finite number of discrete configurations. This formalism 
may be extended to describe systems with continuous 
configuration space, as done in the present note. 



APPENDIX B: EXPLICIT SOLUTION OF THE 
MASTER EQUATION 

We report in this appendix the exact solution, in 
Laplace transform, of the master equation given in 
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Eq. ([T8]), Introducing P(r, s) denned in Eq. Jl9|), the 
master equation becomes: 

1 /"°° Ht' 

S P(r, S )-Po(r) = --P(r, S )+^r) / — P(t',s), 
T Jo T 

(Bl) 

where Pq(t) is the initial condition, Pq{t) = P(t, t = 0). 
Solving for P(r, s), one finds 

P(r, S ) = ^M + ^^), (B2) 

1 + ST 1 + ST 

with tD(s) = J °° dr' t' -1 P(t', s). One first needs to solve 
for uj(s). Dividing Eq. (|B2[) by r and integrating over r, 
one has 



1 + ST 



1 + ST 

(B3) 



Solving for u)(s) yields 

a \Jo 1 + stJ J q 1 + st 
Replacing in Eq. (|B2p . one gets 
tPo(t) 1 t^(t) 



P(t,s) 



1 + ST S 1 + ST 



(B5) 



1 + st' 



1 + st' 



Using the time scale t*(s) introduced in Eq. (|2T|) . the 
solution of the master equation can be written as 



I + ST S T*(s) 1 + ST Jq I + ST' 

The quantity t* (s) may be interpreted as an effective av- 
erage lifetime, in the sense that it is defined as an average 
of t over the distribution iP(t), taking into account a cut- 
off of dynamical origin around a time scale t* « s . In 
the specific case where Po(t) = - 0( r ) I Eq. (|B6|) simplifies 
to Eq. (|20|) after a straightforward calculation. 

Note that it is interesting to introduce a second time 
scale Tq (s) associated to Pq(t) through 



(B7) 



since one can then rewrite Eq. (|B6|) in the more compact 
form 



P(t,s) 



1 



ST*(s) 1 + ST 

x[st:(s)P (t) + (1-st*(s)) V(t)] 



(B8) 



The time scale Tq (s) is the effective average time associ- 
ated with the initial distribution P (t) , taking again into 
account the cut-off time scale t* ~ s _1 . 



In the large time limit, only the asymptotic behavior of 
P(t, s) for s->0 remains relevant. In this limit, the ex- 
pression between square brackets in Eq. (|B8|) reduces to 
-0(t), as st*(s) and stq (s) go to zero (even though t*(s) 
and Tq(s) may diverge when s — » 0), and one recovers 
the same result as Eq. ([20)) : the memory of the initial 
condition is lost. In particular, the scaling form Eq. (j24|) 
and the functional form of the scaling function 4>{x) given 
in Eq. (|2l)]) do not depend on the initial condition. Yet, 



the crossover time t rj s beyond which the scaling form 
Eq. (|24|) becomes valid depends on Po(t). For instance, 
if Po(t) has a peak around a value T p k 3> To, the scaling 
regime is reached only for times t ^> T p k, so that it may 
be necessary to wait for an arbitrarily large time. 

Performing the inverse Laplace transform of Eq. (|24[) 
then leads to the scaling form ([2"5|). where the scaling 
function 4>{u) is given by^ 



<j>(u) 



sin7ra 1 



7rr(a) 



/ dv V 



(B9) 



An asymptotic expansion of Eq. (|B9[) for large and small 
values of u respectively leads to the asymptotic behavior 
of P(t, t) given in Eqs. (|27|) and (|2gj> . 



APPENDIX C: SUGGESTED PROJECT 

The dynamical distribution P(t, t) and the scaling be- 
haviour given in Eq. (|2"rj)) can be obtained using a simple 
numerical simulation of the random process. Drawing a 
random number v from a uniform distribution between 
and 1, compute the lifetime t of the current state as 



T V 



-■L/c 



(the microscopic time tq may be set to 



unity). Check that the probability distribution p(r) sat- 
isfies 



p(r) 



'o 



-l + Q ' 



T > Tq. 



(CI) 



This can be done either numerically, or analytically notic- 
ing that p(r)|dr| = \dv\. Once a state with lifetime t is 
given, the sojourn time t s in this state is chosen ran- 
domly from an exponential distribution of mean r. To 
do this, draw a uniform random number u, < u < 1, 
and compute t s = — Tin it. Check again that t s has the 
desired distribution. Then iterate this computation and 
sum the different sojourn times until the final time t is 
reached (t is an arbitrary fixed time such that t 3> To, 
typically in the range t = 10 3 To to 10 6 t ), and record 
the lifetime t of the state occupied at time t. 

Repeat this procedure a large number of times, and 
compute the histogram of the recorded values of t; nor- 
malizing the histogram yields the distribution P(t, t). 
Compute this distribution for different times, and check 
that for a > 1 , it becomes independent of t in the large 
t limit, whereas for a < 1, the i-dependence remains for 
arbitrarily large t. In the latter case, check the scaling 
form Eq. (|2"rj)) . and the asymptotic behavior of P(t, t) for 
both t <§C t and t ^> t. 
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Alternatively, these simulations may be carried out by the lifetime r. 
considering the momentum p or the energy E instead of 
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